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Abstract 

Background: Parainfluenza virus type 3 (PIV3) was isolated from dairy buffaloes {Bubalus bubalis) naturally affected 
with respiratory and reproductive clinical conditions. 

Results: Examination of nasal and vaginal swabs collected from 12 diseased buffaloes led to the isolation of three 
paramyxovirus isolates from two animals. Antigenic, morphological and biological characteristics of these three 
isolates were essentially similar to those of members of the Paramyxoviridae family. Antigenic analysis by direct 
immunofluorescence and cross neutralization test placed these isolates together with bovine parainfluenza virus 
type 3 (BPIV3). Nucleotide and amino acid phylogenetic analysis of partial matrix gene sequences of the buffalo 
isolates and six field BPIV3 isolates from bovines in Argentina were studied. Buffalo isolates were similar to 
genotype B (BPIV3b) while the six BPIV3 isolates were similar to genotypes A (BPIV3a) and C (BPIV3c). 

Conclusions: This is the first characterization of BPIV3 in water buffalo. 

According to the samples analyzed, in Argentina, the genotype B was found in buffalo and the genotypes A and C 
were found in cattle. 



Background 

Paramyxoviruses are well-known pathogens of the central 
nervous and respiratory systems of many host species. In 
the last few decades, many novel paramyxoviruses have 
emerged causing devastating illnesses in different aquatic 
and terrestrial animals, including in some cases a species 
jump to humans [1]. Parainfluenza virus type 3 (PIV3) is 
an enveloped, single-stranded negative sense RNA virus 
within the Respirovirus genus of the Paramyxoviridae fam- 
ily [2,3]. The Respirovirus genus includes human parainflu- 
enza virus types 1 and 3 (HPIV1 and HPIV3, respectively), 
Sendai virus (murine PIV1) and bovine parainfluenza virus 
type 3 (BPIV3) [2]. In some instances when animals are 
subjected to high stressful conditions, infection with BPIV3 
can contribute to tissue damage and immunosuppression, 
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resulting in severe bronchopneumonia from secondary 
bacterial infections [4]. The resulting disease is part of the 
bovine respiratory disease complex (BRDC) and is consid- 
ered as the most significant illness associated with feedlot 
cattle in the USA [5], and possibly worldwide. Other respira- 
tory viruses such as bovine herpesvirus 1, bovine viral diar- 
rhea virus (BVDV) and bovine respiratory syncytial virus 
(BRSV) have also been associated with BRDC development 
in feedlot cattle [6]. The clinical presentation of bovine 
infections with BPIV3 can vary considerably, ranging from 
asymptomatic infections to severe respiratory illness. In 
most cases where BPIV3 is implicated in disease, usual clin- 
ical signs include coughing, fever and nasal discharge [6]. 

Three genotypes, A (BPIV3a), B (BPIV3b) and C 
(BPIV3c) have been described, based on genetic and phylo- 
genetic analysis [7]. The BPIV3b genotype could hypothet- 
ically be a lineage from a strain that recently crossed from 
another host species into cattle [8]. PIV3 infections were 
found in a wide variety of mammals including cattle, 
humans, sheep [9], goats [10], bison [11], guinea pigs [12], 
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black and white rhinoceros [13], moose [14], bighorn sheep 
[15] and camels [16]. Cross-species infections have been 
reported in numerous instances, including HPIV-3 in 
guinea pigs [12], BPIV3 in a human, BPIV3 in sheep, and 
ovine PIV3 in cattle [17]. In Argentina, serological studies 
conducted in 1980 and 1984 showed a high incidence of 
antibodies against BPIV3 (77%) in cattle from the main 
livestock breeding regions [18,19]. In addition, positive ser- 
ology was reported in domestic and wild South American 
camelids [20-22] and the virus was isolated from cattle and 
sheep [23]. However, little is known about the circulation 
of BPIV3 in cattie and other ruminants in Argentina, 
as well as about the genotypes that are present in this 
Southern American country. 

Production systems have evolved to mixed managements, 
where alternative production species, such as the water buf- 
falo (Bubalus bubalis), coexist in the same habitat as with 
cattle. Water buffalo breeding represents an important eco- 
nomic alternative in Argentina, which allows access to na- 
tional and international markets. This species is susceptible 
to several viral infections [20,24,25] including BPIV3 as 
reported in a very early in 1966 in Egypt [26]. In addition to 
its potential relevance with respect to water buffalo health, 
this finding bears epidemiological significance due to the 
risk of transmission of the virus to cattle. The aim of this 
study was to characterize antigenically and genetically 
BPIV3 isolates from a respiratory and reproductive syn- 
drome outbreak in dairy buffaloes and compare the buffalo 
isolate with contemporaneous BPIV3 circulating in cattle. 

Results 

Morphological and biological properties of buffalo 
viral isolates 

PIV3 was successfully isolated from two out of twelve 
water buffaloes studied, in one case from both nasal (7 N) 
and vaginal (7 V) swabs, and in the other only from the va- 
ginal swab (2 V). 



Isolates produced characteristic cytopathic effects (CPE) 
on MDBK cells, similar to those of the BPIV3 reference 
strain, with many scattered, rounded, refractory cells and 
small syncytia. CPE were first visualized on the second pas- 
sage, after three days post inoculation (DPI), and CPE were 
observed at day 4 or 5 DPI. Both bovine polyclonal anti- 
BPIV3 serum and FITC-anti-BPIV3 IgG reacted with the 
cells infected by these isolates (Figure 1 A). 

Transmission electron microscopy revealed spherical 
to pleomorphic virions, approximately 100-300 nm in 
diameter, which were morphologically indistinguishable 
from paramyxoviruses. Intact virions were enveloped 
and densely packed, as has been described for BPIV3. 
Nucleocapsids were visible and exhibited a typical "her- 
ringbone" pattern (Figure 1 B). 

The three isolates (7 N, 7 V and 2 V) were able to agglu- 
tinate red blood cells (RBCs), with UHA titers of 16, 8 and 
4 respectively. The agglutinated RBCs eluted after 1 h sug- 
gesting neuraminidase function/activity. The HI test 
showed 58% of positive reactions, with titers ranging be- 
tween 1.6 and 2.5 (data not shown). 

Molecular characterization RT-PCR and sequencing 

A fragment of the M gene, consistent with the expected 
size of 328 bp was amplified from the 3 virus isolates (7 N, 
7 V and 2 V) by RT-PCR. As expected, the primer pairs 
did not amplify using genomic material of BRSV as tem- 
plate. Analysis of the alignments of 271 pb of the corre- 
sponding sequences revealed that they are related to 
BPIV3 strain Q5592 from Australia with 94% nucleotide 
identity (Figure 2) and 100% predicted amino acid identity 
(Table 1). 

Data and phylogenetic analysis 

Phylogenetic analysis of the M gene fragment of the three 
Argentinean buffalo isolates (7 N, 7 V and 2 V) together 
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Figure 1 Microscopic examination of PIV3 isolated from buffaloes. A: Positive MDBK cells by IF, as observed by epifluorescence in an 
Olympus BX 40 + H hal microscope (400X magnification). B: Electron micrograph of two purified virions of the 7 N isolate, morphologically 
indistinguishable from paramixoviruses. 
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(See figure on previous page.) 

Figure 2 Nucleotide alignment of M fragment sequences of Argentine buffalo and bovine isolates, other bovine parainfluenza virus 3 
and human parainfluenza virus 3 that circulate in the world. Dots indicate identity with respect to SF strain (AF1 78655). 



with other paramyxoviruses, placed the buffalo isolates in 
the same clade as BPIV3 (Figure 3). 

In particular, the isolates grouped with the BPIV3b geno- 
type represented in the tree by an Australian isolate 
(Q5592). BPIV3b is distinct from the previously character- 
ized members of BPIV3 represented by genotypes A and C. 

In order to determine whether the PIV3 genotype found 
in buffalo is also circulating in cattle, the same gene fragment 
was amplified by RT-PCR using as template BPIV3 RNA iso- 
lated from six bovines of the Buenos Aires region. Three of 
the six isolates show 98% nucleotide identity and 100% 
amino acid identity with BPIV3c (SD0835 strain) reference 
isolate, and the remaining three, 99% nucleotide identity and 
100% amino acid identity with BPIV3a (SF strain) reference 
isolate (Table 1). Phylogenetic characterization showed that 
the bovine isolates belong to the BPIV3a and BPIV3c geno- 
types, which have been reported so far in North America, 
Europe and China (Figure 3). 

Serum neutralization test 

Buffalo sera (n = 24) were analyzed by cross-neutralization, 
Three were negative for both viruses, BPIV3a and b, 9 
showed higher titers against BPIV3b isolated from buffalo, 
4 showed higher titers against BPIV3a and the remaining 
11 had same titers for both viruses. (Table 2). In general, 
sera from buffaloes showed higher antibody titers against 
the isolated buffalo virus than against the reference virus. 

Discussion 

In the present study, we describe the isolation of 
BPIV3 from two water buffaloes with respiratory and 
reproductive signs. The virus was isolated from nasal 
and vaginal swabs from two out of twelve animals 
with reproductive and respiratory signs studied. To 
our knowledge little is known about the pathogenesis 
of this virus in buffalo. Should there be a syndrome 
similar to that present in cattle where the virus is 
transiently excreted [27], the unsuccessful virus isola- 
tion from the then remaining symptomatic animals 
could be related to the time point of sampling, which 
could have taken place at the end of the acute 
infection. 

The IF assay showed strong and specific binding of 
anti-BPIV3 antibodies to cells infected with these 
three positive clinical isolates. Sera from buffaloes 
were able to neutralize the activity of the 7 N, 7 V 
and 2 V isolates and that of the reference bovine 
virus (BPIV3 strain SF), but the neutralization titers 



against the latter were lower than against buffalo iso- 
lates, suggesting antigenic differences. 

The sequence of the M amplified fragment showed 94% 
nucleotide identity between buffalo isolates and BPIV3b 
strain Q5592 from Australia [8]. This is the first report of 
molecular characterization of a BPIV3 isolated from buffa- 
loes. In order to determine homology between PIV3 
isolates from Argentinean bovine and buffaloes, we charac- 
terized retrospectively six isolates collected from bovine 
respiratory outbreaks that occurred in 2009 and 2011. 
These isolates showed high percentages of nucleotide and 
amino acid identity with A and C genotypes (Table 1). 
There are only few reports about BPIV3 circulation in 
Argentina, Epstein et al described the isolation, virological 
and physicochemical characterization of PIV3 from Argen- 
tinean cattle and sheep in 1974 [23]. However, phylogen- 
etic studies of these isolates were not carried out. 
Phylogenetic reconstructions based on the alignment of 
the M-gene nucleotide sequences demonstrated that the 
7 N, 7 V and 2 V isolates clustered in the B genotype 
group. So far, this genotype has only been previously 
reported in Australia [8]. The phylogenetic position of the 
Argentinean buffalo isolates might be correlated with the 
antigenic differences observed by cross-neutralization stud- 
ies, where the reference strain belonging to the BPIV3a 
group was lesser neutralized than the buffalo viral isolates 
by sera of buffaloes. Therefore the serological and molecu- 
lar characterization of buffalo BPVI3 isolates show that 
these viruses are more similar to B than to A (SF reference 
strain) and C BPVI3 genotypes (SD0835 reference strain). 
Interestingly, the genotype B isolated from buffaloes was 
not found among the bovine samples, analyzed in this 
work, that clustered with genotypes A and C. Genotype B 
was until now observed only in Australia with the hypoth- 
esis of a recent crossing from another species into cattle 
[8]. This evidence of genotype B in buffaloes should need 
further investigation to know if this species could be 
one of the host species of this genotype. Therefore ana- 
lysis of a larger number of isolates from cattle and 
buffaloes is needed to determine whether the three 
genotypes circulate in both species, or if the distribution 
observed is representative. Importantly, this is the first 
report of bovine B and C genotype circulation out of 
Australia and China, respectively. The circulation of 
genotype C only in China was attributed to a geographic 
restriction [7]. This hypothesis can be ruled out by our 
results and this observation was most likely due to the 
relative frequencies of these genotypes with respect to 
the samples analyzed. 



Table 1 Comparison of identity percentages in amino acid and nucleotide sequences of the 217-bp region of the M gene among buffalo and bovine isolates 
and reference strains of BPIV3 
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Figure 3 Phylogenetic analysis based on M gene fragment. Phylogenetic analysis was performed using MEGA 4 software with 1,000 
bootstrap replicates. Phylogenetic analysis based on M gene fragment nucleotide sequences compared to other Respirovirus genus members. 
The Argentinean isolates are indicated in bold. GenBank accession numbers are presented in Materials and Methods. 



To our knowledge this is the first report of the iso- 
lation of a parainfluenza virus type 3 from vaginal 
secretions. However, it is important to highlight that 
the presence of the virus in the secretions could be 
due to oro-genital contact between animals. Evidence 
of viral replication in the genital mucosa needs new 
investigations. 



Conclusion 

The isolation and characterization of BPIV3 from 
water buffaloes were reported. We identified the iso- 
lates as bovine parainfluenza virus type 3 genotype B 
that differs from genotypes A and C circulating in 
Argentinean cattle herds reported so far. Further 
studies would be required to determine if A and C 
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genotypes also circulate in buffalo herds, as well as to 
elucidate the role of buffalo as reservoirs, especially 
of the B genotype, and in the transmission of this 
virus to cattle. 

Methods 

Samples 

Nasal and vaginal swabs were obtained from 12 water buf- 
faloes of a dairy farm in the Argentinean province of 
Chaco, during an outbreak of respiratory and reproductive 
disease in July 2009. These animals had white stringy 
mucus in the nose and vulva. These samples were inocu- 
lated on Madin Darby bovine kidney (MDBK) cell cultures. 
The SF strain of BPIV3 was used as reference. Cells were 
daily examined under optical microscope to observe the 
appearance of cytopathic effects. 

Bovine isolates were obtained from cattle with respira- 
tory syndromes, 5 of them in the year 2009 and the 
remaining isolate in 2011. They were collected and ampli- 
fied from Azul laboratories (Azul-Buenos Aires province) 
from outbreaks of Buenos Aires province. 

Immunofluorescence test 

Virus isolates were inoculated into confluent MDBK cells 
in 8-well chambers at a MOI of 1.0, 0.1, and 0.01. For dir- 
ect IF test, Cells with cytopathic effects were fixed with 
99.5% acetone and then incubated with FITC anti-BPVI3 
IgG (Ames). For indirect IF, cells were incubated with anti- 
bovine PIV3 polyclonal antibody (NVSL) diluted 1:100 in 
PBS at 37 °C in a moist environment for 1 h. Then, cells 
were incubated with FITC antibovine IgG (KPL). In both 
cases, slides were layered with buffered glycerin and 
observed under epifluorescence in an Olympus BX 40 + H 
hal microscope. 

Transmission electron microscopy 

Twenty-five ul of viral suspension (16 HAU) were placed 
on a collodion membrane supported on a copper grid, 
negatively stained with phosphotungstic acid, and 
observed with a JEOL 1200 EX II transmission electron 
microscope [28]. 



Hemagglutination assay (HA) and Haemagglutination 
inhibition assay (HIA) 

Serial two-fold dilutions of culture supernatants from 
infected MDBK (50 ul) were titrated against 50 ul Guinea 
pig erythrocytes in V-bottom micro titer plates (Nunc). 
After incubation at 4°C for 4 h, hemagglutination was 
evaluated by the appearance or absence of a red cell 
button. Results were expressed as hemagglutinating 
units/50 ul (HAU/50 ul). The last dilution at which 
hemagglutination was observed was taken as the end- 
point of the hemagglutinating activity. The reciprocal 
of this dilution expresses the number of hemaggluti- 
nating units or virus titer in UHA. 

For the HI assay serial two-fold dilutions of sera from 
buffaloes were allowed to react with a fixed dose of refer- 
ence viral haemagglutinin (8 UHA), followed by the 
addition of guinea pig agglutinable erythrocytes. Titers 
were expressed as the loglO of the reciprocal of the highest 
dilution hemagglutination inhibitory multiplied by a con- 
stant factor of 8. 

RT-PCR and sequencing 

Infected MDBK cells were scraped off the plates and 
homogenized by three cycles of freezing and thawing. 
After an initial centrifugation at 3,000 g for 15 min, 
polyethylene glycol 8000 (Sigma), 10% (w/v, final con- 
centration) was added to the cell lysates, followed by in- 
cubation for 4 h at 4°C. The virus was pelleted at 
12,000 g for 60 min at 4°C. Viral genomic RNA was 
extracted from the virus pellet using RNeasy Mini Kit 
(QIAgen). Oligonucleotide primers for BPIV3 detection 
and identification (Mfwd: 5 'AGTGATCTAGATGATGA 
TCCA 3' nt - 3960 and Mrev: 5'GTTATTGATCCAAT 
TGCTGT -3' nt - 4288) were designed based on a 328 bp 
segment of the consensus BPIV3 Matrix (M) gene. cDNA 
was synthesized using the specific oligonucleotide primer 
(Mfw) and MMLV reverse transcriptase (Promega). 
Subsequently, Go Taq DNA Polymerase (Promega) was 
used to amplify the M gene fragment. 

Mfw and Mrev primers were tested for specificity by RT- 
PCR on BRSV (bovine respiratory syncytial virus) RNA, 
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giving negative results. PCR products were purified using a 
Qiaquick PCR purification kit (Qiagen). Sequencing reac- 
tions were performed twice in both directions with BigDye 
Terminator v3.0 kit (Applied Biosystems, Lennik, Belgium) 
and analyzed with ABI Prism 3730 DNA Analyzer (Ap- 
plied Biosystems). 

Phylogenetic analysis 

M fragment nucleotide and predicted amino acid sequences 
of buffalo isolates were edited and analyzed with BioEdit 
version 7.0.5.3 [29]. Alignments of these sequences together 
with other from representative isolates of previously identi- 
fied BPIV3 and human parainfluenza virus type 3 (HPIV3) 
were created by Clustal W. Nucleotide sequence align- 
ments of the M gene fragment were used for phylogenetic 
analyses using MEGA software, version 4 [30]. Phylogenetic 
relationships were investigated by bootstrap analysis (1000 
replicates) using the Neighbor joining method with the 
Kimura two-parameter model for nucleotide data analysis. 
GenBank accession numbers of the sequences used in the 
phylogenetic tree are as follows: BPIV3 reference strains: 
BPIV3 strain Kansas/15626/84, D84095; BPIV3 strain, NC 
002161, AF178654; SF strain AF178655; Q5592, EU277658; 
SD0835, HQ530153; BPI3 JCU strain, EF 108221; BP4158 
stain EF 108222; BP4159 strain, EF 108223; BP6121 strain, 
EF 108224; BP6128 strain, EF 108225; BP7393, EU266069; 
ISU92, EU439428; Texas-81, EU439429 and HPIV3 refer- 
ence strains: HPIV strains, M16458, M16569; HPIV3-Gpv 
strain, NC 001796 and 14702 strain, EU424062. Sendai 
virus (AB005795) as outgroup. 

Serum neutralization test 

Cross-neutralization tests were performed in MDBK 96- 
well plates, by overnight incubation of 200 TCID 50 on 50 
ul of 7 N buffalo virus isolates or BPIV3 reference strain 
(SF) with 50 ul heat-inactivated two-fold diluted sera. 
The plate was incubated for 5 days at 37 °C, under 5% 
C02, and then screened for the presence or absence of 
cytopathic effects under optical microscope (400x mag- 
nification) to determine the virus-neutralization titers. 
Neutralization titers were calculated as the decimal 
logarithm of the reciprocal of the last serum dilution 
with clear cytopathic effects. 
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